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The 1:1:2 resonant elastic pendulum is a simple classical system that displays the phenomenon 
known as Hamiltonian monodromy. With suitable initial conditions, the system oscillates between 
nearly pure springing and nearly pure elliptical-swinging motions, with sequential major axes dis- 
playing a stepwise precession. The physical consequence of monodromy is that this stepwise preces- 
sion is given by a smooth but multivalued function of the constants of motion. We experimentally 
explore this multivalued behavior. To our knowledge, this is the first experimental demonstration 
of classical monodromy. 
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After more than three hundred years since the formula- 
tion of Newton's laws of motion, one would expect that 
a system as simple as a mass on a spring would have 
been fully understood for some time. In fact, an in-depth 
investigation of even a subset of its possible dynamics 
produces a number of surprises. Chief among these is 
a phenomena known as Hamiltonian monodromy, which 
was introduced by Duistermaat in 1980 as a topological 
obstruction to the existence of global action-angle vari- 
ables [1|. In the resonant elastic pendulum, monodromy 
has easily observable physical consequences. Specifically, 
the observed stepwise precession of the elliptical swinging 
major axis is given by a smooth, but multivalued function 
of the constants of motion. This functional form results 
in loops of values of the constants of motion having dif- 
fering overall behavior, depending on the loop's topology. 

If monodromy were limited to the resonant elastic- 
pendulum system as a special case, it would be consid- 
ered just an esoteric detail. However, it has been shown 
theoretically to exist in many other common and rela- 
tively simple systems, including the spherical pendulum, 
the Lagrange top, and the Kirchoff top .1., ^2, J?, 4]. Most 
intriguing are the quantum mechanical implications in 
atomic and molecular systems. When a classical system 
exhibits monodromy, the energy eigenstates of the corre- 
sponding quantum system can not be mapped onto a sim- 
ple lattice labeled by integer quantum numbers. Defects 
in the lattice of eigenstates are the striking signature of 
monodromy in the global structure of a quantum spec- 
trum. In addition to the static effect that monodromy 
has on global quantum numbers, dynamical consequences 
have also recently been predicted [3, Q . Important quan- 
tum systems have been shown theoretically to have mon- 
odrorny including elhpsoidal biUiards 0, trapped Bose 
gases , the H2 molecular ion , the hydrogen atom in 
combined electric and magnetic fields [10|, [ill, 1 1211 . dipo- 
lar symmetric-top molecules in electric fields [l3|, and 
the ro- vibrational spectra quasi-linear molecules such as 
As we discuss below, monodromy can 




FIG. 1: Diagram of an elastic pendulum in a gravitational 
field of acceleration g. One type of motion is swinging with 
fixed energy in the z-dimension. This motion projected onto 
the XY plane is an ellipse. 



occur near relative equilibria. There is thus the intrigu- 
ing possibility that the singular behavior of monodromy 
may be a common feature of dynamics near chemical iso- 
merization thresholds [l3]- 

A quantum analog of the resonant elastic pendulum 
under consideration here (FiglT]) is the Fermi resonance 
in the CO2 molecule, whose monodromatic features have 
been thoroughly investigated theoretically 
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spite the large number of systems in which monodromy 
is theoretically predicted, there have been no previous 
classical experiments and only a single quantum experi- 
ment fist of which we are aware. In developing a more 
heuristic understanding of monodromy in quantum sys- 
tems, it is useful to have a classical example to guide one's 
intuition. Thus we designed our experiment on a read- 
ily realized classical system in which the consequences of 
monodromy are relatively easy to observe. 

Hamiltonian monodromy is a property of certain in- 
tegrable systems. For a more complete introduction 
than appears here, see Refs. (lil [20 ]. For concrete- 
ness, consider an integrable conservative classical system 
with N degrees of freedom described by the Hamiltonian 
H{qi,Pi) with generalized coordinates {qi} and conju- 
gate momenta {pi}, with i ~ 1, . . . ,N. The {qi,Pi} form 
the coordinates of the 2N dimensional phase space Pq^ ■ 
Such systems contain a set of N independent constants of 
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motion {Fk{qi,Pi)} that have mutually vanishing Poisson 
Brackets, {Fm, Fn}pB = for all m,n. The indepen- 
dence of the Ffe refers to the linear independence of the 
differentials dF^ almost everywhere in Pqp- These types 
of systems allow the construction of local action-angle 
variables [2l[. These variables are a set of coordinates 
(action-angles) and conjugate momenta (actions) that 
obey particularly simple equations of motion. Specifi- 
cally, the actions are all conserved quantities whose con- 
jugate coordinates are 27r periodic and evolve linearly 
in time with constant frequencies. Ratios of these fre- 
quencies, known as rotation numbers, convey information 
about the local choice of action-angle variables. If a ro- 
tation number can not be defined globally by a smooth 
function, then neither can the action-angle variables used 
to calculate it. 

The actual calculation of the set of action-angle vari- 
ables is highly nontrivial for all but the simplest inte- 
grable systems. When monodromy is present, this con- 
struction is even more difficult as the action-angle vari- 
ables are not defined globally. Instead of construct- 
ing these variables directly, the question of whether or 
not a system contains monodromy is typically answered 
through a topological analysis of the system's so-called 
energy-momentum map (EM). The EM takes points from 
Pgp and maps them to the N-dimensional space made 
up of the constants of motion, denoted by . In gen- 
eral, Cp contains two types of points, regular values and 
critical values. Regular (critical) values correspond to 
motions where the differentials dFk are linearly indepen- 
dent (dependent). Examples of motions corresponding 
to critical values are equilibria, with no motion at all, 
and relative equilibria, where motion occurs but the con- 
stants of motion are interdependent. For the elastic pen- 
dulum, relative equilibria include pure springing and cir- 
cular swinging at a fixed height. As a general rule, the 
physically accessible domain of Cp is bounded by critical 
values as the constants of motion are necessarily inter- 
dependent at extreme types of motion. Monodromy can 
occur when critical values are imbedded inside the physi- 
cally accessible domain of Cp . For a concrete example of 
a region of Cp both bounded by and containing imbed- 
ded critical values, see Fig. [H^c). Imbedded critical val- 
ues allow for the existence of loops of regular values that 
cannot be continuously shrunken to a point. After go- 
ing once around one of these singularity-enclosing loops 
of initial conditions, the final action-angle variables will 
differ from the initial ones. Determination of whether or 
not a system contains monodromy then reduces to find- 
ing the constants of motion and analyzing the locations 
and types of critical values for the EM. We now summa- 
rize this calculation for the resonant elastic pendulum, 
the details of which are in Ref. [l^l . 

The three-dimensional (3D) elastic pendulum consists 
of a bob of mass m attached by a spring of constant k to 
a pivot point and placed in a vertical gravitational field 



with acceleration g, as shown in Fig. [TJ Arbitrary initial 
conditions can result in either regular or chaotic motions 
[2^ . We are interested in regular motions resulting from 
small displacements from the equilibrium-hanging posi- 
tion. The two relevant modes of oscillation inherent to 
the system are swinging (pendular) motions with angular 
frequency ujp = y^g/l, and springing motions with angu- 
lar frequency Ug — y^k/m, where I is the equilibrium- 
hanging length. These two frequencies are coupled be- 
cause I depends on k via I — Iq + where Iq is the 
unstretched length of the spring. Particularly interesting 
dynamics occur if the pendulum's swinging and springing 
frequencies are in a ratio 1:2, in which case energy can 
transfer efficiently between the two modes. This trans- 



fer is an example of parametric resonance [2J]. Specif- 
ically, nearly pure springing motion evolves into nearly 
pure swinging motion. The projection of this near-pure 
swinging motion onto the XY plane is an ellipse with 
near-constant major-axis orientation. The plane defined 
by this major axis and z is called the swing plane. Energy 
then transfers back into the springing mode, and this cy- 
cle repeats. The most striking aspect of this motion is 
that the major axis of the projected ellipse precesses in a 
stepwise manner between successive near-pure swinging 
motions (Fig. [3|). Furthermore, the step size between 
the swing planes is constant and a function of the ini- 
tial conditions. Theoretical treatments of this behavior 
have been carried out by Refs. [2^ H^. The main 
result in Ref. [22| is that the stepwise precession A(3 of 



the swing plane orientation /J is a rotation number of the 
integrable approximation. It is given by a smooth but 
multivalued function as a direct physical consequence of 
the existence of monodromy. This function is given in 
terms of the constants of motion, which we now define. 

The construction of the set of constants of motion be- 
gins with an examination of the Hamiltonian. With the 
origin at the pivot point, the full Hamiltonian for the 
1:1:2 resonant elastic pendulum is given by 
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where the unit scaling m — g = I ^ 1 has been used, 
J/ = {uJs/(-iJp) = {kl/mgy/^ = 2 is the ratio of springing 
to swinging frequencies, r — + + and the pi 
are the i*^ dimensional momenta. H is invariant under 
a rotation around the z-axis. Therefore the angular mo- 
mentum Lz — xpy — ypx is conserved. After an expan- 
sion to cubic order about the free-hanging equilibrium 
position and an averaging over the fast dynamics, the 
effective Hamiltonian is given by 

H = \{pI + pI + pI + + y'^ + i^^z^) 

+ (- f [i^P- + ypv) P- + +V^)z- {pi + pI) z] ) 
= H2 + 
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FIG. 2: (a) Loops of constants of motion in (A, x) space gen- 
erated from various sets of initial conditions. The singularity 
at the origin corresponds to pure springing motion, (b) Re- 
sulting swing-plane precession angles for loops corresponding 
to those in (a), to lowest order in (A, x)- Position on a loop is 
labeled by the angles a.i , which are defined from the center of 
the i*'' loop with respect to the vertical. Clockwise is positive, 
(c) Critical values of the EM map consist of a conical lemon 
boundary surface and an enclosed thread. 

where fi = ^{v^ ~ 1): H2 (-ffa) is the quadratic (cubic) 
contribution, and z is now measured from the equilibrium 
hanging position. Under this approximation, the system 
is integrable with three constants of motion {i?, iJ2, L^] 
in three degrees of freedom. The critical values of the 
resulting EM : {x,y, z,px,Py,Pz} ^ {H,H2,Lz} form a 
"conical lemon" surface and a thread [Fig. EIc)]. The 
thread corresponds to purely springing motions at var- 
ious total energies. Loops in the 3D range of the EM 
are somewhat difficult to visualize, so a reduction to an 
equivalent two-dimensional (2D) system is carried out. 
Each set of constants of motion is classified first by its 
value for H2, effectively taking a slice out of the conical 
lemon in Fig. [DJc) . A scaling of the remaining two quan- 
tities {H,Lz} is implemented to map the various slices 
onto each other. This mapping allows different sets of 
constants of motion to be compared easily as points in 
the now 2D lemon [Fig. EJa)]. The new dimensionless 
constants of motion are 

H — H2 Lz 

With this scaling, the thread pierces the 2D lemon at 
(A, x) — (0,0), producing a singularity [Figlll^a)]. This 
singularity is responsible for the existence of monodromy. 

The presence of the monodromy-producing singularity 
causes a rotation number of the integrable approximation 
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FIG. 3: Measured mass positions as projected onto the XY 
plane during three successive near-pure swinging motions. 
These data are fit to the expected elliptical functional form. 
The major axes of these ellipses represent the orientation of 
the swing plane /3. In the example data, A/3 = 32°. 

to be multivalued. This rotation number corresponds to 
the step size A/3 of the stepwise-precessing swing plane 
during a full cycle of swinging, springing, and back to 
swinging. Thus the physical consequences of monodromy 
are easily observed. To first order, A/3 takes the form 

Af3^arg{x + iX). (3) 

The arg function extracts the argument (phase angle) of 
a complex number and can be made single-valued but 
discontinuous through a branch-cut or can be viewed as 
multivalued and smooth. The presence of this multival- 
ued rotation number proves that monodromy exists in 
this system. Explicit calculations for the various loops 
of constants of motion appear in Fig. ^ Curves in Fig. 
^h) are labeled on the right with their corresponding 
loop numbers in (a). Loops not enclosing the singularity 
return smoothly back to their initial values. Loop 5 is 
different in that it does not. 

The experimental goal is to measure A/3 at positions 
along loops similar to those in Fig. [^a) by varying initial 
conditions. Successful experimental measurement of A/3 
relies on relatively pure and long-lived swinging motion. 
The former is necessary to determine when the motion is 
to be classified as purely swinging, purely springing, or 
in transition. The latter allows for one or more complete 
ellipses to be traced out for each purely swinging motion, 
which facilitates a determination of the swing-plane ori- 
entation. The purity and lifetime of the swinging motion 
depends both on initial conditions and on how well the 
2:1 resonance condition is satisfied. Our experimental 
parameters are k — 6.8(2) N/m, m = 0.224(1) kg, and 
unstretched spring length of = 1.00(1) m. These pa- 
rameters yield v = 2.0(1), and pure swinging motions 
that persist for several seconds. Typical energy-damping 
times are on the order of minutes. 

Once motion with well-defined swinging motions can 
be reproducibly created, the next requirement is an ac- 
curate determination of the ball's 3D position as a func- 
tion of time, x(t). These data allow for a determination 
of the swing planes and therefore the stepwise precession 
angle, as well as the calculation of A, x, and the expected 
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FIG. 4: Experimental measurements of the step angle be- 
tween successive swinging motions as loops in (A, x) space 
are mapped out. A loop enclosing the singularity (solid cir- 
cles) does not return to its initial value, thus demonstrating 
monodromy. These data agree with the theoretical predic- 
tion (solid line). A loop [loop 3 in Fig. 2(c)] not enclosing the 
singularity (open circles) is shown for comparison. 



A/?. We capture the motion using two video cameras op- 
erating at 30 frames per second, one (XY camera) shoot- 
ing video from below and the other (XZ camera) from 
the side. To determine the position of the mass, a cir- 
cle is fit to the ball's image for each frame in the video. 
Our circle-fitting scheme relies on the color gradient at 
the edge of the ball's image and therefore necessitates a 
high-contrast ratio between the ball and its background. 
To achieve the needed contrast, the ball is painted white 
and illuminated by intense halogen lights against a black 
background. From measured x{t), we fit the projected 
elliptical motion of the ball, as viewed by the XY cam- 
era during pure swinging motion, to an ellipse using a 
direct least-squares scheme. The major axis defines the 
swing plane. Fitting of successive pure-swinging motions 
yields the stepwise precession. Our procedure is to mea- 
sure three such swing planes and take the average of the 
two steps. Typical data for a single experimental run are 
shown in Fig. [3l Each point represents the fitted position 
of the pendulum bob for a single frame of the captured 
video. Step angles between swing planes {1, 2} and {2, 3} 
typically differed by 2° to 10°, depending on the ellipse 
eccentricity. Instantaneous positions and velocities were 
measured during the part of a projected ellipse with the 
least curvature to minimize errors due to accelerations. 
These instantaneous conditions were used to calculate 
the constants of motion and the expected step angle. 

Step-precession behavior for two experimental loops is 
shown in Fig. [D As in Fig. [D^b), the plots are denoted 
on the right by their approximate loop labels correspond- 
ing to Fig. [U^a). Since swing planes are only determined 
within an additive constant of 180°, continuity of nearby 
ai measurements is used to determine absolute positions 
of the solid circles. Deviation of the singularity-enclosing 
loop (~ 5, solid points) from the expected straight line 
is due mainly to nonideal measurements of the constants 



of motion. In contrast, the noisy structure in the non- 
singularity-enclosing loop (~ 3, open points) is a combi- 
nation of these nonideal measurements and the difficulty 
of experimentally launching the pendulum bob with pre- 
determined constants of motion. This difficulty causes 
the experimental loops to not be perfect circles, but, cru- 
cially, they are still homotopic to circles. This difficulty is 
not an issue for the singularity-enclosing loop, as it shares 
the same origin as the arg function. The qualitative dif- 
ference between the singularity-enclosing loop ('^ 5) and 
the non-singularity-enclosing loop ('^ 3) is clearly evi- 
dent. We see that the former loop does not come around 
to its initial value upon returning to the initial point, 
while the latter does. 

It is a simple distinction: does or does not the size of 
a stepwise precession advance by 360° as one maps out 
the behavior along a loop through constants-of-motion 
space? Yet, the key to our first-ever experimental study 
of classical monodromy is our ability to observe both 
sorts of behavior in the resonant elastic pendulum. We 
hope that this simple classical example can be part of a 
solid foundation upon which to build the intuition nec- 
essary to understand the subtle, but by no means rare, 
instances in which monodromy profoundly influences the 
quantum spectra of atoms, molecules, and more compli- 
cated objects. 
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